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ABSTRACT 

Interior  Search  Algorithm  (IS  A)  is  a  novel  algorithm 
recently  developed  for  solving  optimization  problems. 
Though,  this  metaheuristic  has  proven  to  be  effective 
for  function  optimization  on  widely  used  benchmark 
functions,  this  algorithm  may  converge  slowly  or  be 
trapped  in  local  minima  for  complex  functions.  In  this 
paper,  a  new  hybrid  approach  of  IS  A  is  presented. 
The  proposed  approach  includes  a  fuzzy  logic  system 
for  improving  convergence  of  the  solutions  in  IS  A  by 
dynamic  parameter  adaptation.  To  validate  the 
performance  of  the  proposed  algorithm,  several 
comparisons  are  presented  on  the  basis  of  benchmark 
functions.  Simulation  results  demonstrated  the 
efficiency  of  the  proposed  fuzzy  IS  A  algorithm. 

Keywords —  interior  search  algorithm;  fuzzy  logic; 
parameter  adaptation;  benchmark  functions 

I.  INTRODUCTION: 

Optimization  techniques  provide  an  important 
paradigm  for  solving  real-world  complex  problems.  In 
an  engineering  environment,  optimization  process 
includes  searching  for  a  vector  that  represents  an 
optimal  solution  relative  to  others  for  a  specific 
objective  function.  Optimization  processes  can  be 
divided  into  two  main  categories:  deterministic  and 
heuristic.  The  former  involves  rigid  movements  which 
are  not  adaptable  paper  is  the  proposed  Fuzzy  ISA 
(FISA)  with  dynamic  adaptation  of  parameters  by 
studying  the  effect  of  iteration  count  on  convergence 
and  accuracy  of  solutions  while  updating  the 
parameter. 


This  paper  is  organized  as  follows;  Section  II  presents 
the  ISA  algorithm,  Section  III  presents  the  proposed 
technique,  Section  IV  presents  the  simulation  results 
and  Section  V  concludes  and  presents  future  research. 

II.  INT  ERIOR  SEARCH  ALGORITHM 

The  main  inspiration  for  the  proposed  method  ISA  is 
the  manner  of  arrangement  of  objects  in  the  room  in 
order  to  improve  the  overall  aesthetics  known  as 
interior  designing.  This  method  is  carried  out  in 
accordance  with  the  requirements  and  limitations  set 
by  the  client  [11], 

Keeping  in  mind  the  aesthetics,  the  designer  tries  to 
replicate  the  best  features  of  the  room  which  can  be 
used  as  an  inspiration  for  optimization.  It  may  be 
noted  that  the  algorithmic  procedure  can  be 
mathematically  modeled  by  the  following  processes: 

According  to  the  convexity,  continuity  and 
differentiability  of  the  problem.  On  the  other  hand, 
heuristic  undertake  stochastic  movements  which  are 
flexible,  derivative  free  and  therefore  provide  optimal 
solutions  with  higher  efficiency.  The  main  categories 
include  Evolutionary  algorithms,  Swarm  intelligence 
and  Physics  based  algorithms  [1],  Recently,  heuristic 
methods  inspired  by  human  behavior  have  been 
extensively  studied  and  they  provide  a  convenient 
mechanism  for  dealing  with  different  optimization 
problems.  Some  of  the  popular  algorithms  include 
Teacher  Teaming  Based  Optimization  (TFBO)  [2], 
Tabu  Search  (TS)  [3],  Imperialistic  Competitive 
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Algorithm  (ICA  )  [4],  Harmony  Search  (HS)  [5], 
Mine  Blast  Algorithm  [6],  Seeker  Optimization 
Algorithm  (SOA)  [7],  Social-Based  Algorithm  (SBA) 
[8]  and  Group  Search  Optimizer  (GSO)  [9]. 

Advantages  of  these  algorithms  range  from  low 
number  of  tuning  parameters  [2]  [3],  less  number  of 
operators  than  Evolutionary  techniques  (elitism, 
mutation  etc.)  to  low  dependence  on  few  individual 
solutions  (global  best,  hierarchical  structure  etc.)  to 
guide  the  search  like  in  Swarm  intelligence. 

The  Interior  Search  Algorithm  (ISA)  [10]  is  a  human 
inspired  evolutionary  computing  based  algorithm, 
recently  proposed  by  Gandomi  in  2014.  The  main 
contribution  of  this  Mirror  work. 

In  this  algorithm,  solution  set  is  divided  into  two 
groups,  wherein,  the  composition  group,  the 
composition  of  solutions  is  modified  to  find  better 
solutions  and  the  mirror  group,  where  mirrors  are 
placed  between  the  fittest  element  and  other  elements 
to  find  better  views.  This  division  is  controlled  by  a 
parameter  ‘a’  which  varies  from  a  range  of  0  to  1. 
Since,  ‘a’  is  the  only  parameter  of  the  algorithm,  it  is 
important  for  the  distribution  between  the 
diversification  and  intensification  search  carried  out 
by  the  algorithm. 

A.  Composition  design 

Composition  or  interior  design  follows  a  structured 
procedure  that  includes  a  coordinated  and  systematic 
process.  The  interior  design  usually  starts  from  the 
walls  (bounds)  towards  the  center.  Therefore,  the 
space  becomes  more  and  more  limited  for  designing 
the  composition  of  the  elements  (solutions).  Since  the 
design  space  becomes  limited  at  every  stage  of  the 
procedure,  hence  only  a  more  beautiful  view  (fitter 
solution)  will  be  obtained  at  each  stage  due  to  the 
limitation  s  imposed  by  the  client. 

B.  Mirror  Work 

In  mirror  work,  different  mirrors  are  used  to 
provide  different  views  in  order  to  create  a  more 
decorative  room. 

Hence,  mirrors  are  placed  next  to  beautiful  elements 
(fitter  will  work  as  a  local  search  or  intensification 


otherwise  as  a  solutions)  to  emphasize  their 
effectiveness  (fitness).  This  leads  global  search  or 
diversification  which  can  be  observed  from  to  better 
search  around  the  fittest  element(s)  as  placement  of 
mirrorsnear  global  best(s)  can  lead  to  better  views 
(fitter  As  shown  in  this  section,  ISA  possesses  a  single 
parameter,  solutions).  This  can  be  modeled  by  placing 
a  mirror  randomly  a  which  makes  the 
algorithm  highly  suitable  and  easily  between  each 
element  and  the  global  best  (fittest  solution). 

Therefore,  the  position  of  a  mirror  for  th  element  at 
th  iteration  can  be  formulated  as, 

,  =  3  x  -1  +  (1  -  3  )x  (2) 

where  3  is  a  random  number  between  0  and  1 .  Since 
the  virtual  image  location  depends  on  location  of  the 
mirror,  it  can  be  formulated  as 

=  2x  -  -1  (3) 

Furthermore,  the  position  of  the  global  best  is 
changed  using  a  random  walk  as  local  search.  This 
can  be  formulated  as  , 

-1  +  X  (4) 

where  is  a  scale  factor  depending  on  search  space  and 
is  a  normally  distributed  vector  of  random  numbers. 
The  updated  elements  (global  or  otherwise)  and 
images  (virtual  elements)  are  retained  only  if  the 
fitness  is  improved  compared  to  earlier  position. 


Fig.l  Composition  stages  for  2D  Himmelblau  function 
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Fig.2  Mirror  work  search  part  of  algorithm 


For  an  image  placed  near  the  global  best,  the  mirror 
search  will  work  as  a  local  search  or  intensification 
otherwise  as  a  global  search  or  diversification 

Fig.  2.  The  pseudo  code  of  ISA  is  shown  in  Fig  3. 

III.  DYNAMIC  ADAPTATION  OF 
PARAMETER  ALPHA 

Fuzzy  logic  proposed  by  Prof.  Lotfi  Zadeh  in  1965  is 
based  on  fuzzy  set  theory.  Fuzzy  logic  helps  in 
modelling  knowledge  through  the  use  of  if-then  fuzzy 
rules. 

[iiitialize  the  solutions  randomly  .t,(i  =  1 .2, .  .n) 

Initialize  parameter  a. 
while  (j<  Max  number  of  iterations} 
for  i  =  1  to  n 

Calculate  fitness  and  find  bestxi; 
if  xT  ,  Update  xA  using  Eq  (4) 
else  if  rand  (0,1)  <  a. 

Update  jrj [  usingEq  {2}  &  Eq  (3) 

else 

Update  using  Eq  (1) 

Calculate  fitness  of  nr[T  and  if  better  than 
jsrj ”  1 ,  rep-lac  e  ~ 1  with  jrj 

endfor 
Update  a 
j  =  j+t 
endwhile 


Through  the  use  of  linguistic  information,  fuzzy  set 
theory  provides  a  convenient  method  for  numerical 
computations  by  the  use  of  linguistic  labels  given  by 
membership  functions  [12].  Therefore,  the  proposed 
technique  makes  use  of  a  fuzzy  inference  system 
(FIS)  to  update  the  value  of  a.  In  order  to  properly 
assess  the  runs  of  an  algorithm,  the  iterations  and 
diversity  of  population  need  to  be  taken  into 
consideration.  The  main  metrics  about  ISA  that  are 
used  are  iteration  and  diversity.  Therefore,  the 
aforementioned  metrics  are  used  as  input  variables  for 
the  fuzzy  logic  system  to  adaptively  change  a. 

For  this  purpose,  iteration  is  defined  using  a 
representation  of  percentage  of  populatic$tf{#<e&rf:B£  observt 
be  formulated  as , 

where  is  the  total  number  of  iterations  set  for  ISA  and 
is  the  number  of  elapsed  iterations.  This  metric  is 
used  since  it  allows  us  to  control  the  behavior  of  ISA 
to  perform  a  global  search  in  earlier  iterations 
followed  by  local  search  in  later  iterations.  The  range 
of  input  variable  Iteration  is  from  0  to  1  since  it  can 
only  be  0  -  100%  as  observed  from  Eq.  (4). 
Therefore,  while  starting  the  algorithm,  the  iteration 
will  be  close  to  0  and  can  be  considered  “low”,  and 
when  completed  it  will  be  close  to  1  and  can  be 
considered  “high”.  The  granularity  of  this  metric  onto 
three  triangular  functions  can  be  appreciated  from 
Fig.  4. 

The  next  metric  to  be  considered  is  Diversity  which  is 
a  measure  of  dispersion  of  the  solutions  .  This  can  be 
formulated  as, 

1 

(())=  -Z  X(  -l))2  (5) 

r1  =i 

Where  is  the  current  iteration,  is  the  number  of 
solutions,  is  the  total  number  of  dimensions,  is  the 
size  of  the  population,  is  the  dimension  of  the  solution 
and  the  dimension  of  the  current  best  solution  of  the 
population.  This  metric  is  considered  since  it  can 
control  the  adaptability  of  the  algorithm  to  a  specific 
problem,  with  respect  to  the  diversity  of  the  solutions 
being  obtained.  It  can  lead  to  emphasis  of  exploitation 
when  solutions  are  closer  together  due  to  diversity 
being  “high”  and  exploration  when  solutions  are 
separated  from  each  other  due  to  diversity  being 
“low”.  As  can  be  seen  from  Eq.  (5),  the  range  of 
diversity  will  range  from  a  value  of  0  for  “low”  to  a 
value  of  1  for  “high”. 
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The  type  and  number  of  membership  functions  are 
based  on  experimental  work  carried  out  in  [13].  It 
may  be  noted  that  the  different  type  of  membership 
functions  do  not  show  much  difference  in  results. 

The  output  variable  parameter  a  is  reported  to  have  a 
value  between  0  and  1 .  This  behaviour  was  granulated 
using  five  membership  functions  which  resulted  in 
five  levels  of  “Low”,  “MediumLow”,  “Medium”, 
“MediumHigh”  and  “High”  ,  as  shown  in  Fig.  6.  The 
ranges  that  are  mentioned  for  ISA  that  is  available  in 
the  existing  literature, 


[0.1,0.3] 
[0.1, 0.9] 


Fig.4  Variable  Iteration  as  input  to  FIS 


Fig.5  Variable  Diversity  as  input  to  FIS 


In  the  proposed  technique,  a  is  linearly  increasing  in 
earlier  stages  of  the  optimization  process,  since  a 
lower  a  will  provide  more  exploration  compared  to  a 
higher  a  at  later  stages  of  iterations  which  will  lead  to 
a  higher  exploitation.  This  type  of  division  between 
diversification  and  intensification  is  essentially 
critical  for  any  algorithm’s  performance.  The  range 


[0. 1,0.3]  has  been  claimed  in  [10]  to  be  superior  to 
other  ranges,  hence  this  setting  will  be  followed  for 
conventional  ISA.  The  range  for  FISA  has  been  set  as 
[0. 1,0.9]  for  careful  adaptation  of  parameter  a 
according  to  the  contours  of  the  op  timization 
problem. 


Fig.6  Variable  Alpha  as  output  to  FIS 


Rule  set  used  in  FIS 

1.  If  (Iteration  is  Low)  and  (Diversity  is  Low), 
then  alpha  is  Low. 

2.  If  (Iteration  is  Low)  and  (Diversity  is 
Medium),  then  alpha  is  MediumLow. 

3.  If  (Iteration  is  Low)  and  (Diversity  is  High), 
then  alpha  is  Medium. 

4.  If  (Iteration  is  Medium)  and  (Diversity  is 
Low),  then  alpha  is  MediumLow. 

5.  If  (Iteration  is  Medium)  and  (Diversity  is 
Medium),  then  alpha  is  Medium. 

6.  If  (Iteration  is  Medium)  and  (Diversity  is 
High),  then  alpha  is  MediumHigh. 

7.  If  (Iteration  is  High)  and  (Diversity  is  Low), 
then  alpha  is  Medium. 

8.  If  (Iteration  is  High)  and  (Diversity  is 
Medium),  then  alpha  is  MediumHigh. 

9.  If  (Iteration  is  High)  and  (Diversity  is  High), 
then  alpha  is  High. 
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It  may  be  noted  that  the  five  triangular  membership 
functions  in  the  output  is  due  to  the  fact  that  with  two 
inputs  along  with  their  three  triangular  membership 
functions  results  in  nine  rules  and  operated  on  by 
“and”  operator.  The  rule  set  considered  in  this 
proposed  method  is  shown  in  Fig.  7 

Furthermore,  these  five  membership  functions  along 
with  the  rule  set  represent  the  ideal  behavior  that 
needs  to  be  embedded  in  the  ISA  search  procedure. 
This  can  be  seen  from  the  fact  that  ISA  should  have 
high  exploration  in  initial  iterations,  high  exploitation 
in  final  iterations,  high  exploration  with  low  diversity 
and  high  exploitation  with  high  diversity  of  solutions. 
Therefore,  ISA  performs  a  global  search  initially 
when  the  iterations  are  low  followed  by  local  search 
in  best  area  found  by  algorithm.  The  effect  of 
diversity  allows  ISA  to  escape  local  minima  since  if 
all  solutions  are  in  a  very  small  area  and  not  in  final 
iterations,  then  these  solutions  will  be  forced  to 
separate  by  encouraging  exploration  among  solutions 
due  to  fuzzy  logic. 

It  may  be  noted  that  this  particular  integration  of  FIS 
is  similar  to  the  tuning  of  parameters  as  done  by  an 
experienced  user  of  any  algorithm.  Hence,  FIS  can  be 
effectively  used  to  model  complex  landscapes  of 
optimization  problems. 

The  final  fuzzy  system  is  shown  in  Fig.  8  and  it  can 
be  seen  that  the  FIS  used  is  Mamdani  type  along  with 
the  required  rule  set. 


Fig.8  Fuzzy  Inference  System  for  ISA 


IV.  SIMULAT  IONS  AND  EXPERIMENTS 

This  section  presents  the  used  benchmark  functions 
followed  by  the  obtained  results. 

A.  Benchmark  functions 


The  benchmark  functions  considered  in  this  paper  are 
explained  in  [1]  [10]  [14]  [15]  which  the  reader  can 
refer  to  for  a  brief  introduction.  A  total  of  ten 
benchmark  functions  were  used  for  the  simulations 
consisting  of  six  unimodal  and  four  multimodal 
functions. 

Table  1  presents  all  the  used  benchmark  functions 
along  with  ranges  considered  and  Fig  9  shows  the 
plots  of  these  functions  in  three  dimensions. 

In  all  the  performed  experiments,  the  parameters  used 
for  ISA,  a  is  [0. 1,0.3]  as  in  [10]  and  for  FISA,  a  is 
dynamically  adapted  using  the  fuzzy  inference  system 
discussed  in  the  previous  sections.  The  population 
settings  for  both  were  kept  as  50  for  unimodal 
functions  and  100  for  multimodal  functions.  The  total 
number  of  iterations  is  set  as  1000  for  unimodal 
functions  and  2000  for  multimodal  functions  . 

TABLE  I.  BENCHMARK  FUNCTIONS 


TABLE  I. 

Benchmark  Functions 

Name 

Function 

Range 

Sphere 

C )  =S 2 

[-100,100] 

Rosenbrock 

()=£[*«>£ 

+  (2-1)2] 

[-30,30] 

Quartic 

0  =  Z‘+  [0,i] 

[-1.28,1.28] 

Zackhrov 

+  (£0.5  ) 

[-5,10] 

Step 

0=2(1 +0.4j)2 

[-100,100] 

Matyas 

0=0.26(2  +*)-0*2 

[-10,10] 

Schwefel 

— 

[-500,500] 

Rastrigin 

[-5.12,5,12] 

Griewank 

() 

—  - 

[-600,600] 

Ackley 

[-32,32] 

Since  ISA  is  a  stochastic  process,  a  single  run  is  not 
sufficient  to  verify  the  algorithm’s  robustness  and 
accuracy;  therefore  100  trial  runs  are  conducted  and 
the  results  are  presented  in  the  following  sections. 
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B.  Simulation  Results 

The  simulation  results  are  presented  in  this  section 
wherein  the  best,  worst  and  the  mean  fitness  of  100 
trials  is  presented  along  with  the  standard  deviation 
obtained  by  each  algorithm. 

TABLE  II.  RESULTS  OF  FISA 


Function 

Unimodal  &  Multimodal  Functions  Results 

Criteria 

ISA 

FISA 

Best 

1.7841E-23 

3.6872E-30 

Worst 

2.3818 

1.1584E-20 

FI 

Mean 

0.0795 

4.5839E-22 

Std. 

0.4348 

2.1142E-21 

Best 

6.808  IE-04 

2.3051E-07 

Worst 

3.9867 

4.0119 

F2 

Mean 

0.30231 

0.1681 

Std. 

1.0052 

0.7273 

Best 

2.3339E-04 

5.7577E-04 

Worst 

0.0054 

0.0037 

F3 

Mean 

0.0012 

0.0017 

Std. 

9.2995E-04 

8.0444E-04 

Best 

1.9372E-42 

3.7835E-53 

Worst 

4.5988E-36 

1.0413E-41 

F4 

Mean 

4.2079E-37 

8.0443E-43 

Std. 

1.0688E-36 

2.5161E-42 

Best 

5.9845E-24 

0 

Worst 

0.1281 

5.6633E-22 

F5 

Mean 

4.2796E-04 

1.9967E-23 

Std. 

0.0023 

1.0328E-22 

Best 

7.1402E-182 

1.1362E-318 

Worst 

8.1 179E-173 

1.9457E-304 

F6 

Mean 

3.6344E-174 

7.0168E-306 

Std. 

0 

0 

Best 

-7.1876E+03 

-7.6666E+03 

Worst 

-4.3118E+03 

-4.6486E+03 

F7 

Mean 

-5.8133E+03 

-5.9116E+03 

Std. 

6.9845E+03 

6.8444E+03 

Best 

4.9747 

3.9798 

Worst 

27.8588 

26.8638 

F8 

Mean 

14.990 

13.2329 

Std. 

7.6933 

5.4755 

Best 

0.0344 

0.0295 

Worst 

0.4999 

0.3887 

F9 

Mean 

0.1555 

0.1428 

Std. 

0.0964 

0.0815 

Best 

4.4408E-15 

8.8817E-16 

Worst 

4.4408E-15 

4.4408E-15 

F10 

Mean 

4.4408E-15 

4.3224E-15 

Std. 

0 

6.4863E-16 

As  can  be  seen  from  above  table,  FISA  is  able  to 
efficiently  converge  in  the  case  of  unimodal  functions. 
FISA  achieves  a  much  lower  best,  mean,  worst  and 
standard  deviation  ,  especially  for  FI,  F2,  F4  and  F5. 
It  may  be  noted  that  FISA  is  able  to  converge  to  an 
optimum  of  0  in  the  case  of  F5  which  highlights  the 
high  exploration  as  well  as  exploitation  undertaken  by 
FISA  relative  to  ISA.  In  the  case  of  multimodal 
functions,  FISA  always  performs  better  than  ISA  by  a 
small  margin.  This  is  due  to  the  complex  nature  of  the 
search  space  of  multimodal  functions.  A  notable  point 
may  be  the  performance  of  FISA  on  the  benchmark 
function  F10  wherein  FISA  is  able  to  escape  local 
minima  to  converge  to  a  lower  best  and  overall  mean 
compared  to  ISA’s  performance  on  that  complex 
benchmark  function. 
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Fig.  9  Plots  of  benchmark  functions 


V.  CONCLUSION 

This  paper  presents  an  application  of  fuzzy  logic  to 
improve  the  performance  of  metaheuristic  ISA  by 
dynamically  adapting  the  parameter  of  ISA,  a.  The 
simulation  results  demonstrate  the  efficiency  of  FISA 
in  dealing  with  complex  benchmark  functions  and 
showcase  the  overall  superiority  of  proposed 
technique  of  FISA  over  ISA.  It  may  be  noted  that 
addition  of  FIS  to  ISA  does  not  incur  an  increase  in  s 
pace  and  time  complexity.  Also,  the  use  of  FIS  to 
adapt  a  does  not  lead  to  an  increase  in  the  number  of 
function  evaluations  as  can  be  seen  due  to  addition  of 
operators  such  as  mutation,  crossover,  opposition- 
based  learning  etc.  For  future  research  work,  the 
proposed  technique  FISA  can  be  applied  to  real  world 
applications  like  spring  design,  cantilever  beam 
optimization  etc.  Moreover,  fuzzy  type  2  can  also  be 
explored  in  adapting  the  parameter  a. 
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